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I. 


INTRODUCTION 


It has long been known that the introduction of swirl to a 
flow of reacting fluids can enhance the combustion process. 
Schwartz's experiment, Reference ( 1 ), demonstrated that 
rotating a propane-air mixture in an annular combustion 
chamber results in appreciable shortening of t he flame lengths 
increased flame divergence, improved sta bility chara cteristics 
and delayed blow-off. 


Swirl-can combustors presently under development. Reference 
(2), have also utilized fluid rotation as well as special 
geometric configurations to achieve greater combustion sta- 
bility. 

Detailed measurements of the properties of a turbulent jet 
diffusion flame with swirl have been given by Chervinsky, 
Reference (3). Here too, a shortening of the stabilization 
and flame lengths with increased swirl is observed. An 
accompanying analysis includes a quantitative evaluation of 
the turbulent eddy viscosity and Prandtl number from the ex- 
perimental data for three different degrees of swirl. 

Although improved character!' sties of mixing and combustion 
with swirl have been observed, the basic causes of this 
phenomenon are not well understood. A purpose of this re- 



port is to present a model that accounts for some of the 
effects of swirl as a consequence of the rotational insta 
bility created by its introduction into the flow. 


Mixing length concepts are utilized to demonstrate that the 
addition of a swirl component to the flow enhances the tur- 
bulent transport of momentum. The improved transport is 
evidenced by a modifying factor that increases the eddy vis- 
cosity as a function of swirl and density. 


The model is compared with the three experiments of Reference 
(3) concerning the combustion of a free jet in a quiescent 
medium. It is shown that the formulation 


u a [1 + 90/Np. •) p .u .r . 
T R° J J J 


( 1 ) 


is in good agreement with that data over the range of swirl 
presented . 


In order to extend these results to the case of free jet com- 
bustion in a uniformly moving medium, a computer program was 
developed to solve the governing fluid dynamic equations with 
a flame sheet model incorporated to describe the combustion 
process. After demonstrating that these results collapse to 
the quiescent medium case for low edge velocities, it is 
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found that the eddy viscosity for the jet in a uniformly 
moving medium is modified in a manner similar to that for 
the quiescent medium. 

Since the eddy viscosity model given here implies rapid mix- 
ing when swirl is introduced to the injected fuel, applica- 
tions to confined supersonic combustion are apparent. To 
quantitatively evaluate the effects of swirl on supersonic 
combustion a computer program was developed to calculate the 
finite rate combustion of hydrogen injected from a ring jet 
into an air stream in an annular channel. 

Although mixing is even more rapid for a heavier fuel, hydro- 
gen was chosen because its chemistry is well known. Even for 
the hydrogen case the results of this study show that the 
introduction of swirl to the fuel jet can reduce the com- 
bustor length by one half and result in more uniform burning 
across the width of the annulus. 

It is planned to extend the current results and apply them to 
a variety of new combustor schemes. 
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TURBULENT MIXING MODEL 


II . 


A physical description of the instability of a flow field is 
often given as follows. If a fluid particle is removed from 
its initial position to some other one in the flow field, 
then the background field is stable if a restoring force is 
exerted on the fluid particle. The background field is un- 
stable if a force tends to move the particle further away 
from its origin and neutrally stable if no net force is ex- 
erted on the particle. 


For an axisymmetric flow field with rotational motion this 
concept is easily put in terms of the parameters of that 
field. As shown in Figure (1), a fluid particle originat- 
ing from radius r^ with density p^, swirl velocity w^ is 
considered to have been displaced to a position at r 2 . If 
it is assumed that the process is adiabatic and that the 
radial pressure differences play an insignificant role in 
determining local density differences, then the density of 
the particle remains unchanged at its new position. Since 
there is no torque exerted on the particle its angular mo- 
mentum is constant so that the swirl velocity at the new 
position must be given by 


w 

1 


r i w 


1 W 1 
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FIGURE 1. FLUID PARTICLE DISPLACEMENT IN SWIRL FLOW FIELD 
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II 


The pressure gradient required to balance the centrifugal 
force of the particle at r£ is given by the equation of 
radial equilibrium 


i 2 2 

( dp \ * = Pl( w l) = p 1 r 1 w 1 

\ dr / r 2 r 3 


( 2 ) 


The existing flow field, however, has impressed a pressure 
gradient on the particle given by 



(3) 


Thus, in order for a net restoring force to be exerted on 
the particle, and stability of the rotational flow field, 




p w2r2 > p w 2 r 2 
2 2 2 111 


(4) 


This condition, for constant density, was first given by 
Rayleigh, Reference (4), whereas the above physical descrip- 
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tion has been attributed in Reference (5) to von Karman, 
Reference (6). A mathematical basis for these results has 
been provided in Reference (7) by Synge who arrived at the 
criterion 


d_ 

dr 


(pr 2 w 2 ) =>0 


Stabl e 
Neutral 

Unstable 


(5) 


A motion that is rotationally unstable can result in a state 
of greater turbulence exhibiting increased transport proper- 
ties. In order to determine the effects of rotational in- 
stability on the eddy kinematic viscosity, some essential 
features of the Prandtl mixing length theory i.n a non- 
rotating field are reviewed from Reference (8). 

The eddy kinematic viscosity, e, is. defined by 


e 


dll = 
dr 


I I 

U V 


( 6 ) 


in the case where transport of axial velocity in the normal 
(radial) direction is being considered. If it is assumed 
that the fluctuating components of velocity indicated by the 
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primed quantities, are of the same order, then 


e M a iu'Mv'l a |u ' | | u ' | . (7) 

dr 


As a result of Prandtl's mixing length hypothesis 



(8) 


when SL is the Prandtl mixing length. 

Combining Equations (7) and (8) the Prandtl formulation for 
eddy viscosity is obtained 


e a 


du 
dr • 


(9) 


The eddy kinematic velocity may then be written in terms of 
a turbulent energy transfer as 


e2 a £2 u 1 ^ 


( 10 ) 
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This formulation proves useful for developing arguments on 
the effects of rotational instability on the turbulent vis- 
cosity. 

In a rotationally unstable flow the centrifugal force per 
unit mass impressed upon an eddy is 


(pw 2 ) ' . 
pr 

If it is assumed that this eddy travels a distance i' before 
it arrives at the reference radius then the work done will 
be proportional to 

&'(pw 2 ) . 
pr 

Utilizing mixing length concepts, this expression may be 
arranged in terms of the gradients of the mean background 
flow so that the energy transported due to the rotational 
instability is 


[£ 2 w' 2] a S, 2 ■— ■ (-P- w - ) 
rot. pr 

i n s t . 


sA 9pw 2 . 
pr 3r 


(ID 
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An argument based upon average velocities of mixing (turbu- 
lent fluctuations) can also be given; In principle, it is 
identical to the previous argument but here the details will 
be more apparent. 

Suppose we look at the radius r 0 of the turbulent flow field 
and find some fluid particles moving outward and some inward 
in an eddy-like motion as shown in Figure (2). 

The driving forces on two oppositely moving fluid particles 
should be proportional to the difference between the pres- 
sure and centrifugal forces. Thus, 

.2 2 2 2 
F a P l W l~P.O, W . o) ? a (p 2 w g-P 0 w o) - 
r o 

Assuming the same proportionality, then from Newton's 
second law an expression for the average relative velocity 
between the two particles may be given by 

d ( v j - v 2 ) 3V ^pw* A 2 pw 

dt dt v r 3r a r oP;l " r Q p 2 
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Integrating over some appropriate mixing length, H, and 
replacing the relative velocity by the turbulent fluctuation 
then 


v' 


2 


2 


a 



r o p l 


A£pw* 

r o p 2 


dr. 


To first order this may be written as 


.2 

v 1 


a 


2 2 
(pjWj-^Wg) 


P r 
o o 


& 2 3pw 2 
p o r o 3r 


which is in agreement with the previous result given by Equa 
t i o n (11). 

These results for the effect of rotational instability on 
turbulent energy transport may now be incorporated into the 
eddy kinematic viscosity model. Adding the energy transport 
from rotation effects from Equation (11) to the transport 
from usual shear effects from Equation (8) results in 


( % z u' Z \ +(*. 2 w' 2 

) a I 4 

du 

2 

+ C_ dpw 2 

\ J \ 

'shear 

rot . 

dr 

pr ar 


xm 
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Where C is some constant representing the relative magnitudes 
of the two effects. This constant must be determined by 
experiment. 


An equation for the eddy kinematic viscosity is found from 
Equations (10) and (12) 


£ 


a 





(13) 


The quantity ^ may 


be written 


as 


dpw 2 

"3r 


1 d 
= 7 ? W 


(pr 2 w 2 ) 



(14) 


Since the unstable region requires that (pr^w^)< 0, then 

it follows from Equation (14) that (P w ^)< 0. 

dr 

For this reason the absolute sign appears about this quantity 
in Equation (13). It should be noted that for rotationally 
stable flows there is the possibility that the turbulent 
transport properties of the fluid field can be reduced from 
those of a comparable non-rotating fluid. 
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For free shear turbulence it" has been suggested. Reference 
(8), that overall differences within the flow are more 
significant than local gradients. If this approach is used 
here then Equation (13) can be expressed as 


e 



L o “] l ' Z 

+ C £ 7 |A(pw 2 ) I 


where L is some reference length pertinent to the given 
probl em. 


Defining the turbulent viscosity in terms of the eddy 
kinematic viscosity, and slight rearrangement, yields 


liy = pe = K 


CL | A(pw 2 ) 1 


"I 1/2 


pr 


(Au)' 


PL | Au | . 


(15) 


The form given here is similar to that usually given for 
turbulent free shear flows with the exception that a correc 
tion factor, 


1 + CL. I A(pw 2 ) | 
P r (Au) 2 


1/2 


» 


which includes effects of rotational instability modifies 
the proportionality constant, K. 
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Ill . 


FREE JET COMBUSTION 


(a) Quiescent Medium In this section the eddy vis- 
cosity analysis that has been formulated earlier will be applied 
to the case of swirling free jet combustion in a quiescent medium. 
The consistency of that formulation with available experimental 
data is to be examined. 

It has been shown, Reference (9), that several regimes exist for 
a turbulent swirling jet exhausting into a quiescent medium. As 
the degree of swirl is increased, the adverse pressure gradient 
in the area near the jet is strengthened due to the rapid decay 
of angular momentum. At large degrees of swirl the pressure 
gradient can become large enough to cause a region of reversed 
flow 'near the jet exit. Even before this condition is met the 
gradient can be strong enough to induce a maximum velocity at 
radial positions away from the axis. 

The current work is concerned only with that regime of swirl 
such that within a small fraction of a jet radius downstream 
the maximum velocity appears at the axis. The objective here 
is to predict the effect of this low to moderate degree of swirl 
on the turbulent transport (eddy viscosity) within the jet com- 
bustion region. The only data available for such a correlation 
is that of Reference (3). 
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Observations, References (3) and (10), support the view that in 
cases such as this the swirl velocity is very quickly ( ~ 1 5 radii) 
reduced to one-tenth its original maximum quantity. Yet, the 
effects of swirl are seen to persist to distances of 100 diameters 
apparently through a lasting influence on the transport properties 
of the mixture. For this reason, the corrective factor for eddy 
viscosity of Equation (15) will be based upon initial jet condi- 
tions. Thus, for P e ~P., the correcting factor can be expressed 
by the proportionality 


-1 _ MT 

Re p j u j r j 


1 + C 




(16) 


Here, the jet radius has been used as the length parameter, while 
uniform axial velocity and solid body rotation have been assumed 
for the jet so that Apw 2 = (pw 2 ) evaluated at the jet outer edge 
If the swirl distribution is not solid body then (pw 2 ) max should 
be used. The quantity u/w is often times referred to as the 
Rossby number. For these assumptions and near constant density 
in the jet the Rossby number is related to the degree of swirl 
by 


2 S Nr 0 = 1. 


(17) 
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TABLE I 


CHARACTERIST ICS OF FREE JET OF REFERENCE (3) 


r j = 0 . 082 ft, T j = 300° K , <^- = 5.05, P e = 2116 psf, T e = 300°K 

N Ro 

4.55 

3.02 

1.907 

S 

0.232 

0.160 

0.112 

u j max( ft/sec ) 

249.4 

267.4 

337.9 

w j max( ft/sec ) 

54.8 

88.6 

177.2 

*-s exp/ r j 

11.0 

6.2 

4.0 

*-f calc 7r j 

61 . 

38. 

33. 

L f exp 7 r j 

69. 

51.8 

34. 

2 SN Ro i 

1.02 

0.97 

0.885 


Some of the data of Reference (3) is reproduced in Table I. 
Results have been obtained for only three values of Rossby 
number. The relationship given in Equation (17) is not sat- 
isfied by the case of largest swirl (N Ro =1.907) probably 
indicating that the jet swirl did not conform to solid body 
rotation. The intermediate swirl case (N Rq = 3.02) shows the 
largest discrepancy between calculated and measured flame 


19 












(b) Uniformly Moving Medium In order to obtain 
another viewpoint on the free jet problem, and to gain some 
insight into the appropriate viscosity model for a jet exhaust- 
ing into a moving external stream, a computer program was de- 
veloped to solve the governing equations of the free jet system. 
In this system the exterior region is assumed to have a constant 
velocity. As this velocity is made small, the quiescent jet 
solution should tend to be recovered. 


After applying the boundary layer approximation, the governing 
equations are: 

Conservation of Global Mass 


Tx (pur) + TF (pvr) = 0 


Conservation of Elemental Mass 


pu £_ a k + pv a k = JL 


Conservation of Axial Momentum 


3u 3 u 3n 1^ 

pu 3x + pv 3r " “ 3x + r 


Conservation of Radial Momentum 



k = l ,m 


(19) 


( 20 ) 


( 21 ) 



r 


_ iR 

" 3r 


( 22 ) 
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Conservation of Tangential Momentum 


pu (rw) + pv L (rw) = ± |y 




( r if 


(rw) - 2 rw 


Conservation of Energy 


(23) 


9H 9H 1 9 

pu 3x + P v 9r " r 9r 




1 9H 


N p r 9 r 


F 


+ 1 - 


Np r 


3(u Z /2+w 2 /2 ) 1 v 3_ ] 

^ + J h k N Pr ( N Le ’ 1) 9r a k J 


(24) 


The chemistry of the problem is taken to be satisfied by a 
flame sheet model for the single step reaction of Reference (3) 


| C 2 H 6 + 5C 3 H 8 -+ l C 4 H 10 +5302+200 N 2 
32 C0 2 + 42 H 2 0+ 200 N 2 . 


(25) 


The basic assumption of this model is that at every point in 
the flow field the reaction is completed. Thus, in regions 
where the equivalence ratio, cf> , is greater than unity only 
products and excess fuel are found. Similarly, wherever the 
equivalence ratio is less than unity only products and excess 
air exist. The elemental mass fractions uniquely determine 
the species mass fractions. The required specific enthalpies 
are taken as linear functions of the temperature as indicated 
in Table II . 
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TABLE II 

LINEAR ENTHALPY FUNCTION OF TEMPERATURE 


h k = A k + C p (T-300) 
k 


h k = A k +C p ( T- 2000 ) 
k 


Species T«,(°K) 


C 2 H 6 


c 3 h 8 

c 4 h 10 


(cal/gm) (cal/gm°K) 


(cal/gm) (cal/gm°K) 

1074.0 1.1100 

996.0 1.0980 

1000.0 0.9000 








-2140.0 

0.2020 

| -1640 .,8 

0.32782 

-3205.0 

0.4440 

-2245.8 

0.67856 

0.41 

1 

0.21950 

442.2 i 

j 

0.28216 

0.46 

0.2490 

1 

480.0 

0.3072 


Application of the von Mises transformation to the governing 
equations arranges them into a form more amenable to numerical 
solution. The coordinate transformation is defined by 


pur = i|/ pvr = - t|> 


= 2 f pur 1 dr 1 „ 






















This satisfies Equation (19) and recasts Equations (20) - (24) 
as 


1 Jhe ^uri^k) „ , 

3x - Ijj 3ljJ\Np r ip dip ) » k_1 » m ( 27 ) 


M - 
ax 



(28) 


w 2 _ u ap 
r 2 “ rp dip 


( 29 ) 


lx (ru) ■ | 1(7 _ P T If (n> ' 2r “) 


(30) 



(31) 
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respecti vely . 

The associated initial and boundary conditions for the free 
jet are 


x = 0: u-uj, £ k =a k . H-H JP w = w max 

J j 

\p = 0: M = iil = w = 0 (32) 

3 tp dip dip 

V -*■ °°- u = u e » “k = “k ’ H = H e* w=0, p = p e 


The above equations, along with the perfect gas law 

p = pRT E a k /M k 
k 

properly define the problem. 

A straight forward explicit numerical technique was used to 
integrate these equations. Suffice it to say that the basic 
details of the numerical method can be found in Reference (13) 
with the following modification. 

The inclusion of tangential and radial momentum equations in 
this analysis requires that the normal gradient of the pressure 
be coupled to the axial variation of the dependent variables. 
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This is accomplished by a simple iteration technique. First, 
the pressure is assumed to remain unchanged and the equations 

are integrated one step downstream. New values are obtained 

for the pressure from the latest angular momentum distribution 
by integrating Equation (29) from the edge value to the axis. 

The process is repeated until two successive iterations yield 

axial pressures that agreed to within 0.01 psf. In practice 
no more than three such iterations have ever been required and 
then only for the region very near the initial plane. 

Before applying this analysis to the experiments of Reference 
(3) it is necessary to ascertain the order of the free stream 
velocity that is required to simulate a quiescent exterior 
region. Using an edge velocity of 1% of the jet velocity and 
the two stream eddy viscosity formulation of Kleinstein, Refer- 
ence (11), a computation was made for isothermal air-air mixing. 
Figure (4) demonstrates that the results, are in excellent agree- 
ment with the axis velocity decay profiles of References (8) 
and (11) indicating that a free stream velocity of the order 
of 1% Uj is sufficient to accurately compute axis velocity 
decay . 

This analysis has been used to compare the limited data of 
Reference (3). The eddy viscosity is assumed to be of the 
form given by Kleinstein, neglecting the potential core and 
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with a constant to be determined by matching the experimental 
data 


= K r l/2 l p a u a"Pe u el 

A parameter of great significance for combustor studies is 
the flame length and so this is the criterion that is chosen 
for matching the computations with k of Equation (34). 

Several runs were made with initial conditions corresponding 
to those of Chervinsky (Table I) until the k's for which the 
flame lengths were in agreement could be determined (Table III) 
The flame length here is defined as the distance at which the 
<j>=l line crosses the jet axis. 


TABLE III 

EDDY VISCOSITY CONSTANT AS FUNCTION 
OF ROSSBY NUMBER 


N Ro 

4.55 

3.02 

1.91 

K 

.0315 

.044 

.069 


If the constant is assumed to have the same type dependence 
upon N Ro as indicated earlier, then using the low and high 
N Ro cases of Table III yields an eddy viscosity, 


( 34 ) 
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I 


/ 2 \ 1/2 

P T - .0138 (l + 88.0 /N Ro ) r 1/2 |Pa u a’Pe u el‘ 

This result is in agreement with the variation given earlier 
in Equation (18). In this example it was found that r^ should 
be based upon velocity and not momentum for best results. 

(i.e. r p is that value of r where u = Js(u a +u e )). 

Additionaly, in the current example values of N|_ e = l, Np r = 0.75 
are similar to those obtained by Kleinstein in distinction 
from the analysis of Chevinsky. In the latter work Prandtl 
number was found to vary by more than two orders of magnitude 
through a downstream distance of 100 jet radii. Figure (5a) 
illustrates that the axial velocity profiles using the cur- 
rent matching techniques are good over a large portion of 
the jet decay. 

The axial temperature distributions obtained. Figure (5b), 
are in fair agreement with those of Reference (3), but it 
must be recalled that a huge variation of Prandtl number was 
inferred from those results whereas the current computations 
have assumed a constant Np r< The axial temperature distribu- 
tions given here have been shifted by the distance correspond- 
ing to the stabilization distance of the flame. Reference (3), 


(35) 
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X/ 


ITY DISTRIBUTION ALONG AXIS 






to give this agreement. The general tendency for the computed 
axial temperature to be larger than that observed may imply 
a breakdown of the flame sheet model in this regime as opposed 
to an increased Prandtl number. 

An interesting observation, in agreement with past analyses 
is that within a short distance (x/r,j~15) the swirl velocity 
has been reduced to 10% of its original maximum value. Since 
the swirl correction to the eddy viscosity has been assumed 
to be a function of the initial conditions and the swirl velo- 
city decays so rapidly, then the usual mixing equations with 
the equivalent viscosity are sufficient to attain the perti- 
nent free jet results. That is, in the region of low to 
moderate swirl it is unnecessary to solve the angular momentum 
and radial equilibrium equations. 

To conclude, it appears that an eddy viscosity correction of 
the order (1 + 90/Np o ^) 2 correlates the few existing free 
mixing with combustion data in terms of axial temperature and 
velocity decay as well as flame length. 
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IV. 


AN EXAMPLE - DUCTED RING JET COMBUSTION 


In this section the preceding results are extended to the 
analysis of the mixing and combustion of hydrogen injected 
from a ring jet into a constant area annular air stream. 

The geometry of this system is given in Figure (6). 

The analysis includes the interaction effects of mixing with 
finite rate chemistry. The introduction of swirl into the 
fuel jet is examined with respect to its effect on the char- 
acter of the combustion process. The example chosen is of 
significance for supersonic combustion where rapid mixing 
between fuel and air is required. 


If the problem is considered to be one of free shear mixing 
in a confined region (negligible wall boundary layer effects), 
then the governing equations are similar to those of the free 
jet discussed previously. In this formulation, however, the 
energy equation is written in terms of the temperature and 
chemical species conservation replaces elemental species con- 
servation in order to include the effects of finite rate 
chemistry in the calculations. Thus, 

Conservation of Species Mass 


pu 





NLe 


*Pr 



k = l ,n 


( 36 ) 
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Conservation of Energy 


C D pu IX + TT d pv XI = u XE + uj 

■" 3x H 3r 3x 
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y T N pr 3 r 




( 37 ) 


The heat conduction term of Equation (37) has been modified 
by a turbulent heat transfer term that results from the radial 
pressure gradient, Reference (14). This term accounts for an 
isentropic expansion and contraction of fluid particles fluctu- 
ating between regions of variable pressure. The term is small 
for the present application but is included for completeness. 


As in the free jet case, the equations are more easily handled 
in von Mises coordinates. The transformation here is 


pur = ip X&, pvr = - \p XX. 

8 r 8 x 


(38) 


where 


ip ^ - ^2 _ 2 \ pur 1 dr 1 . 

1 J 

r i 

Since the coordinate, r, is measured from the origin, the 
transformation compels the condition \p^0 for the inner radius 
stream function. The value of ip ^ is arbitrary and was eventu- 
ally chosen in such a way as to enhance the numerical charac- 
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teristics of the computer program. Application of this trans- 
formation results in Equations (28) - (30) given previously and 


ax 



!!k\ \ , 

N pr dr J pu’ k_1 ’ n 


( 39 ) 


c ll = I 3£ + o ul " 2 
P ax p ax ^2 



r2 



2 


N Le £T j^k_ 

a^ ^ c P k dip 


1 1 a c p r f pur 3 T 

pu k W k h k + <p B\p ^TNp r V \p dip 


(40) 



for Equations (36) and (37), respectively. 

The effects of the walls cannot be directly incorporated into 
the analysis because of the different length scales involved 
between the wall boundary layer and free shear layer. Here the 
latter is the major concern and the effects of the wall are 
accounted for by using the results of boundary layer theory to 
obtain a suitable skin friction coefficient, c^ , based upon the 
absolute velocity. The components of the friction force are 
assumed to be in proportion to the components of velocity. An 
a posteriori check showed that the effects of wall shear are 
negligible in the regime analyzed here. The boundary conditions. 
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then are given by 


At 


ip = ^ = o |u = £f ]fcu_/ i + wi\ 

v .34> * dip 2 y T r ^ J 


3/2 


JUaL = ii iUl .2X _ ^ jL =0, r = r.* 
3ib u 94)* dip 2X 1 

pur^Cp 


0 = 4> t 


da 


dip 


^ = 0 


3 u _ c f i|)U / 
dip 2 y T r{ 


3/2 


1 + 


w 


3w _ w 3_u _9T_ ipw c 

dip = u 34)’ 34^ " pur 2c p 


= 0, r = r. 


x = 0: a k = a k (4)), T = T (40. u = u (ip), w=w (4)) 


(41) 


p (0, 4>i ) = P. 


It is not the purpose of this report to detail the computer 
programs developed to obtain the results presented herein. The 
numerical analysis of this problem is similar to that of the 
free jet and follows the standard techniques for ducted flows, 
Reference (15), with few modifications. 

When no swirl is present the annulus outer wall radius is easily 
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matched to the outer wall edge stream function by iteration on 
the axial pressure gradient. All the properties at the next 
step downstream are calculated based upon an assumed pressure 
at new axial location. If the value of r calculated at i|i e is 
not within 1% of the radial mesh width, then the process is 
repeated until convergence. When swirl is present, the situ- 
ation is complicated by a radial pressure gradient. In this 
case the pressure at the next axial station is assumed only 
at the inner radius point. The pressure gradient at a point 
( ip+Aijj , x ) can be related to the pressure at point (ip,x) by the 
expansion 



x,ip + A\p 



x,ip 


3 2 p 

3x3<|> 


AiJj . 


X,1^ 


From Equation (22) 


fih-\ 

r 

23 ( rw) 

rw 3u 4(rw) 3r 

\ 9x3^ J 

u 

3x 

u 3x r 3x 


so that the radial variation of axial pressure gradient can 
be estimated at each step. The iteration now procedes in the 
same manner as for the zero swirl case. 


The mixing model is somewhat simplified here. The jet and free 
stream are both cons i dered, numeri cal ly, i n the initial plane by 
assuming initial profiles as a function of stream function. 
Neglecting the potential core and using the eddy viscosity model 
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of Equation (15), along with the constants developed in Section 
1 1 1 , y i elds 




0.0138 


1 * 9 0( 1 /P) av g t 



w max 

u max 
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for i- 
dr 
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, d 

for j — 
d r 

( pr 2 w 2 ) > 0. 



(42) 


Here the length scale has been taken as the annulus width and 
the average specific volume, (l/p) av g , is 



The model assumes that the effects of swirl are present only in 
the rotati orally unstable regime. 


Coupling of the chemistry with the fluid dynamics of the system 
is accomplished by the method of Reference (16). This provides 
that the one-dimensional chemistry along a streamline may be 
coupled directly to the two-dimensional mixing equations. 

T he hydrogen-air chemistry that is used here, as well as the 
two-step subdomain technique used for its computation, are both 
given in Reference (16). 
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The flow conditions that are investigated here are 


Uj = 12000 fps 
T 1 = 1111. °K 
Arj = 0.11 ft 


u e = 8700. fps 

T e = 1089. ° K 

p(i(j 1 -,0) = 2116. psf 


r. = 1.5 ft 

r e = 2.1 ft 

SI = 0.0025 
2 


( 43 ) 


In the first example the ring jet is placed near the center 
of the annulus width and two cases are presented, one with no 
swirl and the other with Wj = 3300 fps. 


Temperature profiles for these two cases are given in Figure 
(7a, b). It is seen that the combustion with swirl has taken 
place faster and more uniformly. The temperature profile with 
swirl is almost completely uniform at all but the inner radius 
region after 0.37 ft. Tne inner radius region has characteristics 
similar to those of the zero swirl case because the assumption 
has been made that the flow is stable there. This prevents the 
flame from penetrating the fluid between the injector and the 
inner wall. Figure (7c, d) depicts the hydrogen species (H 2 , H, 

OH, H.,0) distribution at x=0.37ft. for the two cases. In par- 
ticular, the H 2 O and H 2 distributions demonstrate that burning 
is more rapid and uniform for the swirling fuel jet. 
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FIGURE 7c, d. 


HYDROGEN SPECIES MASS FRACTIONS AT x=0.37 FT. 






Further indication of the effects of swirl on combustion in 
the annular channel is given by the inside wall pressure dis- 
tributions, Figure (8$. For the case with swirl, the pressure 
rise in the channel is limited to a small axial region, similar 
to a premised case, indicating that the entire mixture is burned 
relatively uniformly. The zero swirl case has a gradual 
increase in the pressure corresponding to a combustion process 
that is still very much coupled to mixing and where burning is 
taking place for some distance downstream. The distance for 
complete combustion appears to have been reduced in half by 
the addition of swirl to the hydrogen jet. 

The combustion with swirl can be made even more uniform by 
introducing enough swirl to the air between the inner wall 
and fuel jet to make this region rotationally unstable with 
respect to the fuel jet. This has the effect of increasing 
the turbulent transport near the inside wall. 

Another method of obtaining better mixing is to locate the 
fuel jet closer to the inner radius. Since the mixing in 
the outer region is more rapid than that of the inner region, 
proper placement of the jet may optimize combustion charac- 
ter i s i tics. 

A first attempt at this was made by locating the ring jet 
about 0.2 ft. inside the location of the previous case. In 
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FIGURE 8. PRESSURE DISTRIBUTION IN ANNULUS WITH 
SUPERSONIC COMBUSTION 




Figure (9a) it is seen that this placement may be too near 
the inside wall since the flame does not reach the outer 
wall as rapidly as before and the inner wall temperature 
has been raised very little. This is probably caused by too 
rich a fuel mixture at the inner wall. 

The jet was then placed half way between the location of the 
two previous cases. In Figure (9b) it can be seen that, al- 
though the combustion length has not been reduced, the temper- 
ature profile at x=0.83 ft. is now uniform over the entire 
width of the combustor. 

Thus it appears that swirl, along with injector location, can 
grossly effect the combustion process by altering resulting 
temperature distributions and the combustion length. 
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IGURE 9. 






CONCLUSIONS 


V . 


A general formulation of some effects of swirl on turbulent 
mixing has been given.. The basis for the analysis is that 
momentum transport is enhanced by turbulence resulting from 
rotational instability of the fluid field. 

An appropriate form for the turbulent eddy viscosity is ob- 
tained by mixing length type arguments. The result takes 
the form of a corrective factor that is a function of the 
swirl and acts to increase the eddy viscosity. The factor 
is based upon the initial mixing conditions implying that 
the rotational turbulence decays in a manner similar to that 
of free shear turbulence. 

Existing experimental data for free jet combustion have 
adequately matched by using the factor (1 + 90/N^ o ) 2 to 
relate the effects of swirl on eddy viscosity. 

The model has been extended and applied to the supersonic 
combustion of a ring jet of hydrogen injected into a constant 
area annular air stream. The computations demonstrated that 
swlrlfng the flow could 

(1) reduce the burning length by one half 
(i'2) result in more uniform burning across the 
annulus width 
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(3) open the possibility of optimization of the 
combustion characteristics by locating the 
fuel jet between the inner wall and center of 
the annulus width. 
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